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Abstract 



We construct a new class of quasi-exactly solvable many-body Hamilto- 
nians in arbitrary dimensions, whose ground states can have any correlations 
we choose. Some of the known correlations in one dimension and some re- 
cent novel correlations in two and higher dimensions are reproduced as special 
cases. As specific interesting examples, we also write down some new models 
in two and higher dimensions with novel correlations. 
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Exactly solvable models have always attracted a lot of interest in theoret- 
ical physics, because they serve as paradigms for more complicated realistic 
models. Examples of such many-body quantum Hamiltonians with exact 
solutions abound in one dimension - e.g., the Calogero- Sutherland model 
(CSM) and its variants P, § - and they have proved useful in the study 
of diverse physical phenomena, such as quantum spin chains , soliton wave 
prop ogat ion || and random matrices |J. A particular feature of these mod- 
els is that their wave-functions show strong inter-particle correlations. This 
feature is expected to be of interest in higher dimensions too, in the study 
of strongly correlated systems. Hence, the construction of exactly solvable 
models in higher dimensions is of great interest. 

A direct generalisation of the CS correlation built into the wave-function 
through the Jastrow factor leads to the Laughlin (or Jain||) wave-functions 
in two dimensions, which are relevant to the phenomenon of the fractional 
quantum Hall effect. Hamiltonians with ultra-short interactions || for which 
Laughlin wave-functions are exact ground state wave-functions have been 
known for years. We have also been able to construct gauge models [|10| for 
which the Laughlin wave-functions are exact ground states. 

Recently, a new form of the pair correlator was constructed in two 
dimensions |TT|1 , which was anti-symmetric under exchange, but also intro- 
duced zeroes in the wave-function whenever the relative angle between the 
position vectors of the two particles was zero or tt. The Hamiltonian for which 
a wave-function constructed from these correlators is a ground state was also 
constructed and studied. Later, this kind of model was generalised to arbi- 
trary dimensions as well [|L2]| . More recently, Ghosh|l3| has constructed new 
CS type models with only two-body interactions in arbitrary dimensions. 
Although these models, a priori, may appear to have little or no physical in- 
terest at the present moment, there are two reasons for studying such models. 
One reason is the paucity of exactly solvable models in higher dimensions, 
which makes such a study a worthwhile exercise. Secondly, and more inter- 
estingly, it challenges the ingenuity of the readers to look for realistic systems 
in physics or in other cross-disciplinary areas, where such novel correlations 
may be realised. 

With this motivation, in this paper, we construct a whole class of exactly 
solvable models in arbitrary dimensions. The ground state wave-function is 
chosen to be an arbitrary homogeneous function (multiplied by an exponen- 
tial factor) of DN variables (D is the dimension and iV is the number of 



2 



particles), which is symmetric or anti-symmetric under exchange of the N 
particles, depending on whether we are interested in bosonic or fermionic 
wave-functions. Then the Hamiltonian is constructed in terms of the par- 
tial derivatives of the wave-function. Further, by writing the Hamiltonian in 
terms of appropriate creation and annihilation operators, a class of excited 
states and their energies is found. 

By taking specific forms of the wave-function, we show that the CSM and 
its variants in one dimension as well as the more recent models in two and 
arbitrary dimensions are reproduced as special cases. Interestingly, we find 
that in two dimensions, we can construct exactly solvable models for which 
the more natural Jastrow-type correlations are relevant. More specifically, we 
show that our procedure allows for the construction of a model for which the 
unprojected Jain wavefunctions are exact ground states. As a final example, 
we construct a new fermionic Hamiltonian in three dimensions with two and 
three body interactions which can be exactly solved for the ground state and 
and a class of excited states. The ground state wave-function of this model 
has the novel feature that it vanishes not only when two particles coincide, 
but also when their position vectors differ by multiples of a specified lattice 
vector (generalising the result of Murthy et al[ll] in two dimensions). 

Let us consider a system of N particles in D dimensions with position vec- 
tors I"i = (x\, xf, • • • , ?2 — x \i ' ' ' i • E 2 > )> • • • i r N = { x Ni x %i ' ' ' i x n)- 

We now prove a general mathematical theorem which is applicable to this 
system. 

• Theorem I : Let / be a function of ND variables (i.e., f = f{xl), i = 
1, 2, . . . , N, j = 1, 2, . . . , D) that is homogeneous and of degree n. Then 

ip — f v (xl)e 2 ^ i ^ Xi isQ an exact eigenstate of the Hamiltonian 



1 N - 1 N a, N D d 2 f/dx j2 

i i i j J 

+ f ED^) ! (i) 



2 Note that f v — g is also a homogeneous function of its variables. Hence, the most 
general form of the theorem has T) = 1. However, for future convenience in constructing 
explicit examples, we retain rj. 
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when 

gi = r] and g 2 = r](rj - 1) (2) 

and with eigenvalue 

E = ^(2r)n + ND). (3) 
Proof : Considered as a function of a single variable x\, we see that 

&f = — 7^ — ^ + — 7 — ^ 

- 2 V xi^^iP-iP + 4 2 ip. (4) 

Summing over all the iVZ) coordinates and using Euler's theorem for 
homogeneous functions, 

N D df N df 



■ j <).i' i 



we get 



z2wj2 = ^-i)E — p — ^ + ^E — - f — i> 

i,j {J x l ij J ij J 

N 

- 2n7]ip - NDijj + ^ r 2 ip. (6) 

i 

Using this in Eq.([T]), we get 

Z i,j 1 

+ 3 (gi ~ ty) E % V> + |(2^ + ND)^. (7) 

Hence, for the choice of g\ and g% given in Eq.^), ip is an eigenstate of 
the Hamiltonian with the eigenvalue given in Eq.(|3T). 
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The next step is to prove that ip is actually a ground state wave-function. 
For this, we formally introduce creation and annihilation operators as follows 



... df/dx 
= «. -* + «!— 

.+ . . df/dx 

A x = Px + ix - if] — j — (8) 

for each of the coordinates x 3 i . We then find that the Hamiltonian in Eq.(|l|) 
can be written as 

1 N 1 
H=^( A l i A x i ) + ^Vn + ND). (9) 

hi 

ip is annihilated by all the Aj and is therefore the ground state of the Hamil- 

i 

tonian with the energy (2r]n + ND)/2. Note however, that [A X ,A\] for dif- 
ferent coordinates x and y do not commute for non-zero 77 and hence, the 
Hamiltonian is not fully soluble for the excited states. 

We can also look for excited states of the Hamiltonian by making the 
following ansatz[ll] for the excited state wave- function - 

= M4)<K4) (io) 

- i.e., we assume that the excited state wave-function factorises in terms of 
the ground state wave-function ipo and a completely symmetric wave-function 
<p. This form can be justified by looking at the asymptotic properties of the 
Hamiltonian. 

Using the ansatz in Eq.fllOD in the original Hamiltonian, we find that 
4>{x\) satisfies 

i,j UJ ^i i,j UJj % i,j J <yj^i 

Here, we have assumed that 

H$ = E m $ (12) 



following ||11|| . We now specialise to the case of radial excitations. - i.e., we 
look for solutions where is only a function of t = J2fj x i 2 - We then find 
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that the eigenvalue equation for <p reduces to the confluent hypergeometric 
equation given by 

+ {Eo-t)^ + \{E m - £ o )0 = 0. (13) 

The solutions are confluent hypergeometric functions which are normalisable 
provided 

^(E m -E ) = m (14) 

where m is a positive integer. This gives the energy spectrum of radially 
excited states as 

E m = E + 2m. (15) 
Note that we have obtained the results for the ground state wave-function 



and the excited states in analogy with Ref.|TTJ; however, our results are much 
more general and are valid for any homogeneous function /. In fact, for the 

function ip = f v e 2 ^ * to qualify as an iV-particle wave-function, 
f v must be symmetric/antisymmetric under exchange for bosons/fermions. 
This is a requirement that we need to impose on /. Furthermore, there exists 
no restriction as to dimensionality, number of particles, types of correlation, 
or whether there exists two-body, three-body or iV-body interactions in the 
potential. But for useful results, we need to put some constraints on /, such 
that either the type of correlation present in the wave-function or the type 
of interaction present in the potential is specified. In either case, we can 
construct a whole host of models with exact ground states. But before we 
go on to the construction of new models, let us first check that by choosing 
appropriate forms of /, several earlier known models are reproduced as special 
cases within this formalism. 

• (I) One dimensional models (D=l) 

— a) Let us choose / to be of the form 

N 2 

/ = [[(xi — Xj) \Xi — Xj\ a and tp = fe ^ 2 . (16) 

i<j 

This is the so-called Jastrow form and is a highly correlated wave- 
function. It picks up a phase (— 1) A under exchange of two particles 
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and can be chosen to be fermionic or bosonic by choosing A = 1 
or A = respectively. By differentiating /, we can check that 

i f i<j ( X i ~ X j) 2 

and the three-body term vanishes, since (without loss of general- 
ity), we may choose x± < x 2 < ■ ■ ■ < x^ and use the identity 

N 1 

E = 0. (18) 

i^jj^k \ Xi X j)\ X i X k) 

Thus for this value of /, the Hamiltonian in Eq.(|I|) yields the 
Calogero- Sutherland model. 

— b)A simple generalisation is to choose / to be of the form 

JV 

/ = n>< - x j) x \ x i - x j\ a ( i9 ) 

i<j 

and construct the Hamiltonian. (k = 1 is the CSM.) For k = 2, 
the three body-term vanishes [fL3fl . But for other values of k, we 
have both two-body and three-body interactions. 

In principle, we can construct many further models with multi-particle 
correlations, since / is arbitrary. For fermionic wave-functions, / has 
to be antisymmetric under exchange, which can always be achieved by 
choosing / to be of the form of a determinant multiplied by symmetric 
factors. 

(II) Two dimensional models (D=2) 
Here again, there are several possibilities. 

— (a) As an obvious generalisation of the one dimensional Jastrow 
form, we may take 



\7l) ll X «l l X jl 



N 

f = II(N-I X ; 

i<j 

and if) = fe (20) 
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where Xj = (x},xf) and / is antisymmetric under exchange when 
A is an odd integer. The Hamiltonian is easily constructed just 
as in the one dimensional case. The three-body term vanishes by 
the same logic that it vanishes for the one- dimensional case - i.e., 
since |xj| is a number, the phase space can be split into identical 
copies and within each copy, it is easy to show that the three-body 
interaction term vanishes by the identity in Eq . flT8|) . 

(b) Recently, a new type of pair correlator was found in two dimen- 
sions with which a Jastrow type correlator could be constructed. 
This was of the form XiUj — x^y^. Using this, one can construct 
the many-body wave-function 

N 

f = T[( x iyj - xjyi) x \xiyj - X0i\ a (21) 

i<j 

and the corresponding Hamiltonian, which was earlier written 
down and studied in Ref . [[TT|| . 

(c) We may also choose / to be of the form 

N 

/ = n(|x i | fc -|x J f) A .||x i | fc -|x i | fc | a (22) 

In this case, using our general procedure, we can write down the 
explicit models. Once again, since |xj| is a number, it is easy 
to check that the three-body term vanishes for k = 2 (besides 
k — 1, of course) and hence for this case, a CS type model can be 



constructed. This model was recently constructed by Ghosh [[13 
as an example of a CS type model in higher dimensions. 

(Ill) Three-dimensional models (D=3) 

Here, again, some of the lower dimensional models can be obviously 
generalised. 

— (a) A direct generalisation of type (c) models in two dimensions 
yields the wave-function 



^ = J](| Xi | fe - |x j | fe ) A .||x i | fc - \xj\ k \ a e~^^ Xi (23) 



s 



for which Hamiltonians can be constructed. As before, for k = 1 
and k = 2, three-body terms vanish. 

(b) A generalisation of the type (b) models has also been earlier 
considered 0]. The simplest way is to consider that the two vector 
determinant is now replaced by the three vector determinant, i.e., 



Pi, 



i 



Vi 

y.i 







Vi 


Zi 




Xj 


Vj 


z i 






Vk 


z k 



P; 



ijk 



(24) 



(Xi 



Hi) in two dimensions and Xj = (xi,yi,Zi) in three di- 



mensions). / can now be constructed as / 



{Pijk) \Pijk\ 



Some models have been considered in even higher dimensions, but we shall 
not consider them explicitly, since our point is merely to illustrate that all 
models with polynomial correlations can be obtained from our general pro- 
cedure. 

Finally, we construct new models in two and three dimensions using our 
general procedure. 

• (i) D=2 

Choose ip to be the unprojected Jain wave-function given by 



N 



\2m„ 



l * 
2L z ^ 



(25) 



where we have used complex notation and \ n is the Slater determinant 
of n filled Landau levels - i.e., Xn is a determinant involving at most 
n—1 powers of z* in each of its terms. For this Jastrow-Slater form, the 
wave-function does not factorise neatly in the usual Jastrow form. Still 
by using / = H^ f < j(zi — Zj) 2m x n and taking the appropriate derivatives 
as in Eq.|l|, we can construct the Hamiltonian, which, however, has 
many body interactions and is not simple. 
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(ii) D=3 

Choose / to be of the following form 



JV 


1 


x i 


Xj 


X 


1 


x i 


Xj 


/ = n 


1 


Vi 


Vi 


x abs 


1 




Vi 


i<3 


1 




Z 3 




1 




Z 3 



N 



i<3 



and tjj = / e -3Ef(^+w?+^) 



(26) 



A fermionic model can always be constructed by choosing A to be 
an odd integer, since the determinant enforces antisymmetry under 
exchange. However, note that the wave-function vanishes not only 
when the relative angle between the two position vectors is zero or 
7r (as in R,ef.||ll||, but whenever the two position vectors satisfy = 
mYj + k(l, 1, 1) (to, k are integers) due to the property of the determi- 
nant. This wave-function obviously has only two-body correlations and 
can be thought of as yet another generalisation of case (b) in two dimen- 
sions. The Hamiltonian can be constructed from this wave-function by 
our general procedure and is given by 



H 



1 v 
2 



N 



(a + \) x . ( yj - z 3 f + {z,j - Xj ) 2 + [ Xj - yjf 



+ 



2 

(a + A) 



E 

2 iV N 

-E ( E 



x ij 



(Vj -Zj), 2 



N 



i j^i 



+ (E 



\ Z 3 X j)\2 



+ (E^-^) 2 

j+i A « 



(27) 



This Hamiltonian clearly has both two-body and three body interac- 
tions. Its energy is given by E = ((a + X)N(N - 1) + 3N/2). 
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(iii) Generalisations 

Simple generalisations of this model with two-body correlations where 



(28) 



1 


X i 


Xk 




1 X? 


x k 


1 




Vk 




i y? 


y n k 


1 




Zk 




1 zf 


~n 
Z k 



or models with three particle correlations where 



Xi Xj Xfc 

y% Vj Vk 



Z.; 



z k 



Vi 
z r - 



X 3 

Vi 



x k 
Vk 



(29) 



are also possible. 



However, our aim is not to write out an exhaustive set of models. It is 
more to demonstrate that our general procedure allows the construction of 
explicitly soluble Hamiltonians. For any specified correlation of the ground 
state wave-function, we have a systematic procedure to construct the Hamil- 
tonian. And for all these Hamiltonians, we can construct a class of excited 
states based on radial excitations, using the ansatz in Eq. (|10|) . 

To conclude, in this letter, we have constructed a class of non-trivial 
Hamiltonians, which can be solved exactly for the ground state and a class 
of excited states. The ground state can have any correlation we choose. In 
particular, we find that any polynomial correlation including the Jastrow 
and several recent novel correlations can be reproduced as special cases. We 
also construct new exactly solvable fermionic Hamiltonians in two and three 
dimensions as an illustration of our general procedure. 
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